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Abstract. We study travelling-wave solutions for a reaction-diffusion sys- 
tem arising as a model for host-tissue degradation by baeteria. This system 
consists of a parabolic equation coupled with an ordinary differential equa- 
tion. For large values of the 'degradation-rate parameter' solutions are well 
approximated by solutions of a Stefan-like free boundary problem, for which 
travelling-wave solutions can be found explicitly. Our aim is to prove the ex- 
istence of travelling waves for all sufficiently large wave-speeds for the original 
reaction-diffusion system and to determine the minimal speed. We prove that 
for all sufficiently large degradation rates the minimal speed is identical to the 
minimal speed of the limit problem. In particular, in this parameter range, 
nonlinear selection of the minimal speed occurs. 
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1. Introduction 

In this article we continue our mathematical analysis of a model for the degra- 
dation of host tissue by extracellular bacteria. This model was introduced in 
[KKC+03] and consists of a reaction-diffusion equation coupled with an ordinary 
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differential equation. In [HKR07j we proved the existence of solutions to the 
time-dependent problem and the convergence to a limit problem in the 'large- 
degradation-rate' limit. Here we turn to the question of existence and behaviour of 
travelling-wave solutions. 

There is an increasing interest in models which support the understanding of 
bacterial infections, and we refer to [KKC+03 ; and !WKK+04] for further back- 
ground and references. This paper is in effect concerned with the specific issue of 
how rapidly a bacterial infection in, for example, a burn wound may invade the 
underlying tissue (with dire potential consequences for the patient, notably mor- 
tality due to septicemia). For the type of model with which we are concerned 
here, the relevant invasion speed is expected to be governed by the corresponding 
travelling-wave problem. Typically, the smallest possible wave-speed is realized by 
a large class of solutions. Accordingly determining the minimal speed of travelling 
waves becomes a central question (with obvious implications for the amount of time 
available for medical treatment, for instance). 

In a dimensionless form, the model in KKC"'"03| is given by the equations 

dtu = /!^u — u + w — ^ku(l — w) , (1-1) 
dtw = ku{l-w), (1.2) 

where u describes the concentration of degradative enzymes, (1 — w) the volume 
fraction of healthy tissue and 7, k are positive constants. The key parameter here 
is the degradation-rate k, which is very large in practice. Equations (jl.l[) . (|1.2p are 
considered in a time-space cylinder, with the upper half space of as the spatial 
domain. Finally, the system is complemented by initial conditions for u and w, a 
Neumann condition on the lateral boundary for u and a decay condition for u and 
w in the far field. In |HKR07| we gave a precise mathematical formulation and 
proved the existence and uniqueness of solutions to a slightly more general system, 
including the possibility of a diffusion term in (|1.2|) . One noteworthy aspect of (II. 
()1.2|) is the convergence of solutions to the solution of a Stefan-like free boundary 
problem as the degradation rate k tends to infinity. This large-degradation-rate 
limit was identified by a formal asymptotic analysis in [KKC^OS] and was proved 
in |HKR07j . 

Reaction-diffusion systems of the general form 

Vt = DAV + F{V), (1.3) 

where V : (0, T) x R" ^ is vector- valued, F : -> is a given nonlinearity 
and D E R^^^ is a diagonal positive-semi-definite matrix, appear in a lot of 
different scientific areas. One-dimensional travelling waves are solutions on (0, 00) x 
R of the special form 

V{t,x) = V{x-ct), 

where c is called the speed and V the profile of this travelling wave. The question 
of existence and behaviour of travelling waves is of enormous interest in many of 
the applications and pertinent results for the vector case N > 2 remain restricted 
to rather specific systems. 

The system (fLT]) . (fOj) has, as we will see in Remark 1 2.1) one stable equilibrium 
in {u,w) = (1,1) and one unstable equilibrium in (u^w) = (0,0) and therefore 
belongs to the class of monostable systems. Scalar monostable equations 

Vt = Av + f{v) (1.4) 
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where 

/(O) = /(I) = 0, ,f{v) > for < f < 1, /'(O) > > /'(I), 

are well-studied, especially the famous Fisher equation, that is (|1.4p with /(w) = 
v{l — v), introduced in |Fis37| . The rigorous analysis of equations of this type 
also started in the 1930s with the work of Kolmogorov, Petrovskii and Piskunov 
[KPP37| . Under an extra assumption on / they proved the existence of travelling 
waves for all speeds c > cun, where cim can be found explicitly in terms of / by a lin- 
earisation about V = (corresponding to a degenerate node in the travelling-wave 
phase plane). Moreover, they proved that the solutions to ()1.4p with initial data 
decaying sufhciently fast propagate with speed cun- For more general monostable / 
the propagation speed was found to be either equal to or larger than cun and there- 
fore one distinguishes between a linear or nonlinear selection of the propagation 
speed (the terminology pulled and pushed fronts, respectively, having an equivalent 
meaning). Aronson and Weinberger [AW75| (see also Hadeler and Rothe [HR75| 
and [Sto76j for other pioneering work on such matters and [vSOSj for a recent re- 
view) proved that for general monostable /, in both the linear and the nonlinear 
selection case, the propagation speed for solutions with initial data decaying suffi- 
ciently fast is given by the minimal speed of travelling waves. They showed that 
monotonic travelling waves exist for all speeds c > Cmin and none for c < Cmin, 
where Cmin > cun; solutions of (II. 4p with sufficiently rapidly decaying initial data 
propagate with speed Cmin- For the nonlinear selection cases, Cmin > Qin, Rothe 
[Rot81| and Roquejoffre |Roq97| proved that, if the initial data decays sufficiently 
rapidly, the large-time solutions to (|1.4p not only propagate with speed Cmin but 
also approach the profile of a travelling wave with minimal speed. 

Whereas the connection between large-time behaviour and existence of travelling 
waves for monostable equations is satisfactorily resolved, the distinction between 
nonlinear or linear selection is still a challenging question, see for example |vS89j . 
[EvSOO| , |BD94] . Only a few rigorous results for general monostable nonlinearities 
are available. In |LMN04| a variational characterisation of travelling waves and a 
concrete criterion for whether linear or nonlinear selection occurs for a given equa- 
tion was derived. Even fewer general analytical results are available for monostable 
systems. An existence theorem for travelling waves was proved in |VVV94] for 
monotone monostable systems, which are systems of the form (jl.3p in which the 
Jacobian matrix DF has only nonnegative off-diagonal elements. Results on the 
existence of travelling waves and the long time behaviour of (|1.3p for monostable 
gradient systems, that is for n = N and nonlinearities F = with g : R" — > R, 
were given m |Mur04j . To the best of our knowledge there are no more general 
results on the question of whether linear or nonlinear selection will occur. 

In this article we prove that for all 7, fc there exist monotone travelling waves for 
the system (II. 1|) . (|1.2p for all speeds c > Cmin and no speeds c < Cmin- The minimal 
speed Cmin in general depends on the parameters 7, k. 

We prove that for all k > ko, where ko is explicitly given in terms of 7, the 
minimal speed is larger than the value Ciin(fc) obtained from a linearisation at 
the unstable equilibrium. Surprisingly enough, for k > ko the minimal speed of 
travelling waves is identical to the minimal speed of travelling waves for the Stefan- 
like limit problem that was formulated in [IIKR07| . Our analysis is based on two 
main facts. One is the monotone structure of the system, which makes possible the 
use of comparison principles for the parabolic problem. The second is a remarkable 
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reduction in order of the travelling-wave equations that occurs when c is given by the 
minimal speed Coo of travelling waves of the large-degradation limit. We obtain the 
existence of travelling waves with speed Coo for (|l.ip . (|1.2p . Finally, a comparison 
argument allows us to prove that nonlinear selection occurs for sufficiently large 
values of fc, the minimal speed in this regime being identical to Coo- 

This paper is organised as follows. In section [2] we prove the existence of travel- 
ling waves for the reaction-diffusion system (jl.ip . (|1.2p . In section |3] we recall the 
formulation of the large degradation limit problem and consider travelling-wave 
solutions for this problem. In section |3] we return to the reaction-diffusion system 
and investigate the selection of the minimal speed. Section [5] deals with the con- 
vergence of travelling waves for the reaction-diffusion system to travelling waves of 
the Stefan-like free boundary problem as the reaction rate k approaches infinity. 
Finally, we give some conclusions and remarks on open problems in section |6l 

2. Existence of monotone travelling waves 

In this section we prove the existence of monotone travelling waves. First we fix 
some notation and make some remarks. 

Remark 2.1. A system of the general form (|1.3p is called monotone if the off- 
diagonal elements of the Jacobian matrix DF are non-negative and strictly mono- 
tone if they are positive (see |VW94j ) . A system of the form l|1.3p with two 
stationary points is called monostahle if one of the stationary points is stable and 
the other is unstable. 

The system ([LT|) . (fOI) is of the form (fL3| with 

It follows that (|l.ip . (|1.2p is a monotone but not strictly monotone system. Further, 
z = (0, 0) and z — (1,1) are the only stationary points of (|l.ip . (|1.2p and we obtain 

with one positive and one negative eigenvalue, and 

with two negative eigenvalues. Therefore (jl.ip . (|1.2p is a monostable monotone 
system. 

As remarked before, for a one-dimensional travelling wave u^w of (II. ip . ()1.2p 
with speed c, the functions 

M(t,a;) := u{x — ct), w{t,x) := w{x — ct), 

are solutions of (|l.ip . (|1.2p on (0,oo) x M. Therefore u,w have to satisfy the 
travelling-wave equations 

= u" + cu' -u-\-w-'-fku{l-w), (2.1) 

= cw' + ku{l-w), (2.2) 

We restrict our investigations to functions u, w taking values only in [0, 1], which is 
the range of meaningful values in the tissue degradation model. We first summarise 
some basic properties of travelling-wave solutions. 



TRAVELLING WAVES IN A TISSUE-DEGRADATION MODEL 



5 



Lemma 2.2. Assume that u e C^(R), w e C^M) with {)<u,w <l satisfy (f2?T|) . 
()2.2|) wit/i c G R. T/ien f/ie following properties hold. 

(1) Invariance under space shifts. For any xq €z R- the functions defined by x t-^ 
u{x + xq), X I— !■ w{x + a;o) solve (|2.ip . (|2.2p wii/i i/ie same c G M. 

(2) Invariance under inversion. The functions defined by x t-^ u{—x), x w{—x) 
solve (|2.ip . (|2.2p wiit/i c replaced by — c. 

(3) Smoothness of travelling waves. The functions u^w are infinitely differentiate . 

(4) Monotonicity of travelling waves. If c > and u,w are not both constant then 
u,w are strictly monotone decreasing and they approach unity as x ^ —oo and 
zero as X ^ CO. 

(5) Standing waves. For c = either {u, w) = (0, 0) or {u, w) = (1, 1). 

Proof. The statements (1), (2) are immediate from (j2.ip . (|2.2p . Statement (3) 
follows from a bootstrapping argument. 

To prove (4) we observe from p.2p that w' < and w' < if u(l — uj) > 0. Next 
we differentiate (|2.ip and obtain 

= u'" + cu" - {1 + -ik{l - w))u' + {I + -iku)w' . (2.3) 

From this equation we see that u' cannot have a positive local maximum. Since u 
is bounded, u' cannot approach a positive supremum at +oo or — oo and we obtain 
that u' < 0. 

Assume now that u'{xq) — 0. Therefore u'{xq) is a local maximum of u' and we 
deduce from (12. 3p that w' [x^)) > and that, using w' < and (j2.2p . 

k 

u(l — w)(xq) — w'(xq) — 0. 

c 

Let us consider the case u{xq) = (the case that w takes the value one is analogous). 
Then, by w' < and w > 0, we obtain u = on [xo,oo). From (|2.ip it also follows 
that u; = on [xo,oo). Since u,w are smooth we find that u,w solve (12. ip . (|2.2p 
with m(xo) = u'(xo) = w(a:;o) = 0. On the other hand solutions of (|2.ip . (|2.2p 
with these data prescribed at xq arc unique and thus identically equal to zero, in 
contradiction to our assumption of non-constant solutions. Therefore u(l ~ w) > 
and u',uj' < holds. 

Next, since u{x),w{x) are monotone in x and uniformly bounded, their limits as 
X — oo exist. From (|2.2p we deduce that w'{x) has a limit and, from < w < 1, 
that 



w'{x) -> 0, m(x)(1 
From (|2.ip we similarly obtain that 
(it' + cm)' (a;) 

and thus 



— w{x)) -^0 as X ^ — oo. (2-4) 
— > as .T ^ — oo. 



—u{x) + w{x) — > as X ^ — oo. 

By < u,w < I, ([2^ . and m', w' < this gives 

u{x) — + 1, w{x) —^l as x — *■ — oo. 

The proof that u{x), vu{x) approach zero as x — > +oo is similar. 
To prove claim (5) we first deduce from c = and (|2.2p that 

u{l -w) = 0. (2.5) 
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Let US assume that u{xq) > for a point xq G M. Then there exists a maximal 
interval (a, b) such that xq E (a, b) and u > in (a, b). From ()2.ip . ()2.5|) . and c = 
we then obtain that 

u > 0, w = 1, u" < in (a, 6). (2.6) 

If a > —oo we deduce from the regularity assumptions on u, w and u > that 

u(a) = 0, w{a) = 1, u'(a) = 

and finally u"{a) > 0, which yields a contradiction to (|2.ip . This shows that 
a — ~oo; by analogous arguments we obtain that 5 = oo and (a, b) = M. Therefore 

= u" - u + 1 on M. 

and from the boundedness of u it follows that u = 1 ~ w. 

Similarly we prove that w = = u if there exists G M such that w(xq) < 1. □ 

Due to Lemma 12.21 we can restrict our investigations to the following set of 

admissible functions. 

Definition 2.3. Let JC be the set of functions v e C°°(R), with 

< V < I, v' < 0, (2.7) 
v(x) — + 1, as X —oo, (2-8) 
v(x) — > 0, as X ^ +00. (2-9) 

We call (c,u,w) G M x /C x /C satisfying ^2.1\) . i2. 2\) a monotone (decreasing) 
travelling wave for (II. 2|) . 

We observe from ()2.2p that the speed of a monotone (decreasing) travelling wave 
is always positive. 

Tiieorem 2.4 (Existence of travelling waves). For all 7, k there is a positive number 
Cmin < OO; Cmin = Cniin (77^)7 such that there exists a monotone travelling wave 
of (|1.2p with speed c for all c > Cmin and such that there is no monotone 

travelling wave with speed c < Cmin- 

The value Cmin gives thus the minimal speed of travelling waves for (|l.ip . (|1.2p . 
In the remainder of this section we prove Theorem 12.41 

For < e < (27)""'^ we consider the following strictly monotone and strictly 
parabolic approximation of (|1.2p . 

dtUe = dxxUe - Ue + We - iku^i.'^ - Wf.), (2-10) 

dtWe = edxxWe + e{ui. - li^) + ku^{l ~ lii,), (2-11) 
and the corresponding travelling-wave equations 

= + CM^ - Me + We - 7fcUe(l - We), (2.12) 

= ew'^ + cw'^+e{ue~Ws) + kueil-We). (2.13) 

Observe that the ODE system (|2.12p . (|2.13p has the same stationary points as 
the original system fM]). For e = the systems (|2T2l) . (|2T3l) and ([2T|) . 

p.2p coincide. The existence of travelling waves for the auxiliary problem and a 
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variational characterisation of the minimal speed follows from |VVV94| . To state 
their result, we define functionals $i, $1 : ^ x /C ^ M, 

g"{x) - g{x) + ajx) - jkg{x){l - a{x)) 
<Pi[g,a) := sup — , 

xm -g[x) 

ea"{x) + e{g - a){x) + kg{x){l - ajx)) 

$1 e,cr sup — . 

xeR -cr'(x) 

Lemma 2.5 ( [VVV94 ] Theorem 1.4.2). For all c > c^, where Ce > is defined by 



Ce ■■= inf max $f(e,cr), a) , (2.14) 

there exists a monotone travelling wave (c, and is the minimal speed of 

travelling waves for (PTTU| . ^(TU) . □ 

To prove the existence of travelling waves for (jl.ip , (|1.2p we pass in (I2.12p , (|2.13p 
to the limit e — > 0. First we derive bounds for u^, which are uniform in e > 0. 

Lemma 2.6. Consider e > and let [c^u^^w^) he a monotone travelling wave for 
(|2.10p . (|2.1ip . Define fi^iy > to be the positive solutions of 

^? + c^l-l = 0, (2.15) 

z/^ - cj/ - (1 + 7fc) = 0. (2.16) 

Then 

u'e > -Ai(l-Ue), (2.17) 

u'^ > -vu., (2.18) 

Kl < c^+l + 7fc, (2.19) 

< > -^(1-zz;.) (2.20) 
c 

hold. 

Proof. We fix a: G M and define for y e (— oo, x\ 
u{y) := l-(l-«,(x))e^(«-^), 

u;(y) := 1 - (1 - w,{x))e^'^y"'\ where L := -^-tf. 



Then we observe that 



lim u{y) = 1, u{x) = ^^(.t), 
lim w{y) — 1, i()(x) = We{x). 



Moreover 



ii" + cu' — u + Wg, — "fku(l — w^) 
= -(1 - Ue{x)){^? +cfi- l)e''(^-^) - (1 - w,) - 7/cw(l - We) < 

and 

ew" + cw' ~ ew + eu^ + fcue(l — w) 
= -(1 - u;e(a;))(£i^ + cL-{e + ku,))e^'^y-''^ - £(1 - u,){x) < 0. 
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By the maximum principle, which we apply once for the scalar equation (|2.12[) and 
once for the scalar equation ()2.13p . we deduce that 

< M, < w on (— oo,x], 

thus 

> u'{x) = — — u^{x)), 

w'^{x) >w'{x)=-L{l-w,{x)), 

which proves (^17)) . The estimate ([^1^ follows from ((^1^ and (PTf|) . 

In addition, by comparing on [x, oo) with 

:= Ue(x)e-''(^-^) 

we obtain 

u < on[a;,oo), (2-21) 
u\x) < u'^{x), 

which yields (|TT5)) . □ 

By similar arguments one proves that corresponding properties hold for 
travelling- wave solutions of 11111), pr2)l . 

Lemma 2.7. Let (c, m, lu) 6e a monotone travelling wave for p.ip . ()1.2p . Define 
fijV > to be the positive solutions of 

H^ + cfi-l = 0, (2.22) 

z/^ - ci/ - (1 + 7fc) = 0. (2.23) 

Then 

u' > -fi{l-u), (2.24) 

u' > ~uu, (2.25) 

\u"\ < c^i+l+-/k, (2.26) 

ix;' > --{1-w) (2.27) 
c 



We prove now the first statement in Theorem 12.41 
Proposition 2.8. For each c > cq, where 

Co := liminfce, (2.28) 

there exists a monotone travelling wave {c,u,w) for p.ip . p.2p . Moreover, the 
value Co is finite. 

Proof. Assume first that co < oo and fix an arbitrary c > co and a subsequence 
£i ^ (i ^ oo) such that 

Co = lim Ce^. 

i — *oo 

By Lemma 12.51 there exists a sequence of monotone travelling waves (ci,Ue. 
for (|2.10p . (|2.1ip . with Ci = c if c > cq and Ci = Cg. if c = co such that 

= u". + CiU^^ - Ue, + Uie, - 7fcUe,(l - WeJ, (2.29) 

= £jtf" + Ciw'^. + ei{ug- — Wei) + kue-{l — WeJ. (2.30) 
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Since travelling waves are invariant under space shifts, we can assume without loss 
of generality that 

iie,(0) = i foralHeN. (2.31) 

By < u^.^w^. < 1, the monotonicity of w^. and (|2.17p . (|2.19p there exists u G 
Ci'i(R),w e L°°{R) such that for all < a < 1, R>0 

Ue. ^ u in Ci'"(hi?,i?]), (2.32) 
Wei ^ pointwise almost everywhere in R (2.33) 

hold for a subsequence et 0{i — > oo). 

Multiplying ([i?29l) . ([2?30| by a function rj e C;?°(M) and integrating we deduce 

= / ( ~ ^'"ii ^'^VK, ~ ^("si - "^Ei ) - VlkUsi (1 " We , 

^ /(.,'..-e,,V...,K,-.,,..,.,,:-..,). 
Due to (|2.32p . (|2.33p we can pass to the limit ^ in these equations and get 



= 



^ — rj'u' + crju' — rj{u — w) — rj^ku{l — w)^ , 



= / ( ~ c?]'w + kriu{l — w)) . 
Jr 

It follows that {c,u,w) solve (I2.ip . (|2.2p and, by a bootstrapping argument, that 
u,w are smooth. Moreover, (|2.3ip and (|2.32p yield that u{0) = 1/2 and by Lemma 
12.21 we obtain that c > and that {c,u,w) is a monotone travelling wave. Since 
c > Co was arbitrary, the first part of the Proposition is proved. To prove that 
Co < cxD, let g, a be two smooth strictly monotonically decreasing functions with 

g{x) ~ = e^^ for all a; > 1, 

g(^x) — a{x) = 1 — e^^ for all x < —1. 

Then there exists a constant C(k) < oo such that for all < £ < 1 and all |a;| > 1 
g"{x) — g{x) + a{x) — 'ykg{x){l — <j{x)) 



-g'{x) 

ea"{x) + e{g — (t){x) + kg{x){l — (t{x)) 



< Cik), 

< C{k). 



-a'{x) 

We can estimate the same ratios for all < e < 1 and \x\ < 1 by a constant 
depending only on ||£i||c2([-i,i])j l|o'llc2([-i,i]) • By Lemma [^751 and the definition of 
Cs in (|2.14p it follows that 

Ce < C{k) for aU < e < 1. 
In particular, co < C{k) is finite. □ 
We now complete the proof of Theorem [ 



Proposition 2.9. There is no monotone travelling wave for p.ip . (jl.2p with speed 

c < Co. 
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Proof. Assume that {c,u,w) with c < cq, u,w £ K. satisfies (|2.ip . p.2|) . We then 
obtain that 

<PUu,w) = c, 

IP/ N w" (x) + u(x) — w(x) /„„,x 

< c + esup — j/- —. 2.34) 

:reR -w'(a;) 

From (|2.2p and w < 1 we deduce that 

u-w u{l-w) c 

T ^ = 7- (2.35) 

—w —w' k 



,u' 



A differentiation in (|2.2|) yields that 

—cw" — ku'{l — w) — kuw' — —cw' kuw' 

u 

which gives, together with p.25|) . that 

w" —u' k . 
= u < V 2.36) 

—w u c 

with ly as in Using in (PIM)) we obtain that 

^l{u,w) < c + eC, 

where C is independent of e > 0. From Lemma [^31 and the definition of cq, (|2.14p 
and (|2.28p we deduce that 

Co < liminf max(<i>^(u, w), $2('"i w)) = c, 

which is a contradiction to our assumption c < cq. □ 

3. The Stefan-like limit problem 

The reaction-diffusion system (jf .ip , (|1.2p converges to a Stefan-hke free boundary 
problem as k tends to infinity, see |HKR07| . For solutions Uoo , M^oo of this limit 
problem 

= C/oo(l-l^oo), (3.1) 

holds and the spatial domain splits in a region where Uoo = 0, and a region where 
Uoo > and Woo = 1- If we denote their common boundary at time t by T{t) then 
Uoo, Woo satisfy 

dtUoo = AUoo - f/cc + 1 in {Woo = 1}, (3.2) 

jdtWoo = Woo in {Uoo = 0} (3.3) 

and a continuity and jump condition on T(t), 

[Uooit, .)] ^ 0, (3.4) 

-[VC/oo(t, .) • Ht, •)] - l[Woo{t, .Mt, •) • Ht, ■), (3.5) 

where v(t, .) and i/(i, .) are the velocity and the unit normal of the free boundary 
T{t), pointing into {Uoo — 0}, and [.] denotes the jump across the free boundary 
from the region {Uoo > 0} to {Uoo = 0}. 

As for the reaction-diffusion system, travelling- wave solutions Uoo , Woo are given 
by a speed c G M and profile functions Uoo , Woo , 

Uoo{x,t) = Uooix-Ct), Woo{.X,t) = Woo{x-ct). 
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We are interested in monotone travelling waves which connect unity and zero, 

Uoo{x),Wooix) 1 as x ^ —cxD, (3.6) 
Uqo{x),Woc{x) — > as a; — > oo, (3-7) 
Uoo,Wqo are monotonically decreasing. (3-8) 

Due to the shift invariance of travelling waves, the condition p.ip and the mono- 
tonicity of , Woo we can assume that 

Uooix) = for x > 0, (3.9) 

Uooix) > 0, Wooix) — 1 for X <0. (3.10) 

From ([X^ - llXS)) we then obtain 

= + + 1 fora;<0, (3.11) 

= 7cw^ + Woo for a; > 0, (3.12) 

and the continuity and jump condition 

Uoo(O-) = 0, (3.13) 
jc(l-w^{0+)) = -u'jQ-). (3.14) 

Proposition 3.1. For all c> Coc, where 

(3.15) 



v/7(l + 7) 

there exists a unique solution Uoo,Woo of (|3.6p - p.l4p . This solution is given by 

~ ' (3.16) 

f 1 tfx<0, 
Woaix) = < _^ (3.17) 
[Pe ifx>0, 

where 



a = a{c):=-^ + ^- + l, (3.18) 
P = P(c):=l-^. (3.19) 



7C 



For c < Coo , there does not exist any solution. 



Proof. We deduce that (PH]) and (EH), (jXT^ hold if and only if holds, with 

a as in p.lSp . Similarly p.l2p and (|3.7p are satisfied if and only if Wrx, satisfies 
([XT7)1 with /? e R. Since 

1 - i«oo(0+) = 1-13, -u'oo(O-) = a 

we obtain that the jump condition (|3.14p is satisfied if and only if (3 satisfies (|3.19p . 
Finally, by (|3.17p and p.l9p . the solution Wao is nonnegative if and only if c > 

Cno- □ 
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u{x) , 









\wix)J 




^0 


_k 
c 



4. Selection of the minimal speed 

In this section we prove that for sufficiently large values k > ko sl nonlinear 
selection principle determines the minimal speed of travelling waves. The threshold 
ko is obtained explicitly in terms of the given constant 7. First, we have to analyse 
the behaviour of travelling-wave solutions at infinity. 

4.1. Behaviour at infinity. The linear selection principle for the minimal speed 
is based on the analysis of the linearised system at the unstable stationary point 
{u,w) = (0,0) of (|2.ip . (|2.2p . In the next lemma we show that solutions have to 
decay exponentially to zero as x tends to infinity. 

Lemma 4.1. Let {c,u,w) be a monotone travelling wave. Then 

lim —logu{x) ~ lim —\ogw{x) = A, (4.1) 

X^QO X x— >oo X 

where X is a negative root of the cubic equation 

= X^ + c\^ -{l+-fk)X--. (4.2) 

c 

Proof. With the definitions 

f-c 1+jk -1\ 

I , Bo := 
--^ / 

the system (|2.ip . (|2.2p is equivalent to 

V^{x) = AoVoix) + ku{x)w{x)Bo (4.3) 

and the linearized system at (0, 0, 0) is given by 

e{x) = Ao^x). (4.4) 

The eigenvalues of ^0 ai'e the solutions of (|4.2p . Since s 1-^ s'^ + cs^ — (1 + 7/0)5 — -I is 
negative at s = and becomes positive as s — s- +00 there exists a positive eigenvalue 
A+ > 0. The other two solutions of (|4.2p satisfy the equation 

= A2 + (c + A+)A + ^ (4.5) 

cA+ 

which has, depending on the values of c, A:,A+, either two negative roots, one re- 
peated negative root or two complex-conjugate roots with negative real part. Since 
Vo{x) converges to zero as x -\-oo, the curve x t—^ Vo{x) is for sufficiently large 
values of x contained in the stable manifold and converges exponentially to zero; 
see for example jPerOl] Section 2.7. By [CL55| Theorem XIII. 4. 5 there exists a 
solution ^ = (^1,^2,^3) of the linearized system (|4.4|) and S > Q with 

Voix) = C(^)+0(e("-^)^), (4.6) 

ax) = Oien 

as X —I- +00, where a < is the real part of an eigenvalue A of Ao. One checks that 
Ao has no eigenvector with a component equal to zero; therefore 

^^{x) = 0(e"^) fori = 1,2, 3 (4.7) 

holds as a; ^ -l-oo. Let us show that in fact A is real. Assume that X = a + if3 with 
P ^ 0. Then ^(x) describes, as a; ^ -l-oo, a spiral around the origin contained in 
the plane P spanned by the real and imaginary part of an eigenvector of Ao with 
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eigenvalue A. But then, since the difference between ^ and Vq decays exponentially 
faster than Vq has to take values outside the set {z = {zi,Z2,Z3) : 21 < 0, 22 > 
0, 23 > 0}, which is a contradiction to the assumption that (c, u, w) is a monotone 
travelling wave. This shows that A = a < is real. Thus (|4.ip follows from (|4.6p 
and (113. □ 



The equation (14. 2|) connects the speed and the decay rate at +00 of a travelling 
wave. We now further analyze this relation. 

Lemma 4.2. For all k > and all X < there exists a unique value Cfe(A) > 
such that A satisfies (|4.2p for c = Cfc(A). The function A 1— s- Cfe(A) attains a positive 
minimum cun = ciin(fc) at a unique value Xim — Aiin(fc) and cim, Aun are given by 

Alin < 



Moreover 



with 



3Afi„ = -(l + 7/j) + 2V(r+7fcF+3fc, (4.8) 

AiinQin = (l+7fc) - V(l + 7fc)' + 3fc. (4.9) 

c;(A) < /orA<Aun, (4.10) 

c;(A) > /or A > Alin, (4.11) 



Cfc(A) ^ 00 as A — > —CO, 
Cfc(A) — > 00 as A ^ 0. 

Proof. From (|4.2p we obtain that Cfc is given by 



c-.(A) ^ 4(^-(^+^'^)l)+Vi(^-(^+^^)l)' + X^' 
that Cfe is strictly positive and that Ck tends to infinity as A — > — cxd or A — > 0. 
Therefore the positive minimum cun is attained at a value Aun and c'(Aiin) ~ 
holds. By (g^l) this imphes 

= 3AL + 2AiinCiin- (1 + 7^)- (4.12) 



Moreover, by (g^), 

= Ajl„ + cii„Afi„-(l + 7/c)Aii„-— . (4.13) 

Clin 

One checks that (j4.8p . (|4.9p is equivalent to (|4.12p . (|4.13p . In particular, cj, has 
only one zero and we deduce that cJ,(A) < for A < Aun and c^(A) > for 

> A > Alin. □ 

Corollary 4.3. The minimal speed Cminik) of travelling waves satisfies the estimate 

Cmin(fc) > Cli„(fc). (4.14) 



Proof. According to Lemma 14.11 for a monotone travelling wave with speed c, a 
negative root A of (|4.2p exists. On the other hand ciin(fc) gives the minimal value 
of c such that (14. 2p has a negative solution. □ 



We now state the result corresponding to Lemma [4.1l if we consider x approaching 
-00. 
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Lemma 4.4. Let {c,u,w) be a monotone travelling wave and let k>l. Then 

lim - log (l - u(a;)) = fj,, (4-15) 

1 k 

lim - log (l - w(a;)) = -, (4.16) 



-oo X C 



where fi > is given by 



c c 

M = -i, + \-r + ^- (4.17) 




/ -u'{x) \ 




^7 


1 


1-uix) ], 


.iH 







\l-w{x)J 




^0 






2 V 4 

For < k < 1 the limits in (|4.15|) . (|4.16p are either equal to or equal to k/c 
Proof. With 



Bi 



(|2.2p is equivalent to 
Vlix) = AiVi{x) + k{l-u{x)){l^w{x))Bi (4.18) 
and the linearized system at Vi — (0, 0, 0) is given by 

C{x) = AiC(a;). (4.19) 
The matrix Ai has the positive eigenvalues 




c I ^ k 

and one negative eigenvalue 



/^2^-^-y^ + i. 

Since Vi{x) — s- as a; ^ — oo we deduce that Vi{x) is for sufficiently small x 
contained in the unstable manifold of (|4.18p at (0,0,0). Using [CL55j Theorem 
XIII. 4. 5 we obtain the existence of a solution ^ of the linearized system (|4.19p and 
a (5 > with 

Vi{x) = C(a;) + 0(e('3+*)"), (4.20) 
Cix) = ©(e'^^) (4.21) 

as a; — !■ — oo, where (3 — k/c ov [3 — ^. One checks that k/c > fi holds for fc > 1. 
Thus, if /3 — k/c then the trajectory of Vi as a; ^ — oo has to be tangential to 
the eigenspace of Ai corresponding to the eigenvalue k/c. On the other hand this 
eigenspace is spanned by a vector (fc/c, l,a), where a < for k > 1. Therefore 
Vi leaves the region {z — (21,22,^3) : 22,23 > 0}, which is a contradiction to 
u{x), wlx) < 1. This proves that /3 = /i for fc > 1 and implies that the trajectory of 
C{x) for X —00 is tangential to the eigenspace corresponding to the eigenvalue 
/i, which is spanned by the vector {fi, 1, 0). Therefore 

C2{x) = Oief^n 
and imni) follows from (g^Ol)- By (1^ we deduce 
d k 

— \og(l ~ w{x)) = - + 0{e^'') as X — > -00, 
dx c 
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and thus (|4?T6l) holds. □ 

Corollary 4.5. Let k>l. If we consider two different monotone travelling waves 
for ()1.2|) . the one with the lower speed converges faster to unity as x ap- 

proaches — oo. 

Proof. For a monotone travelling wave {c,u,w), by Lemma [4.41 the convergence of 
w, w to unity as a; ^ —oo is exponential with rate ^ for u and rate k/c for w, where 
II is given by (|4.17p . Since 



we deduce that both convergence rates are decreasing with c. □ 



4.2. The reduced system for c — Coo- It was observed in KKC"'"03 that the 
travelling- wave equations (|2.ip . (|2.2[) are for the speed c = Coo defined in p.lSp 
remarkable in being equivalent to a system of two first-order equations. 

Lemma 4.6. Let u £ 1Cq,w G /C. Then {coo,u,w) is a monotone travelling wave 
if and only if u, w satisfy 

u = ^Cao{u~ w), (4.22) 
w' = -fc7(l + 7)coo'«(l - w). (4.23) 

Proof. Multiplying by c^o and using we see that (li?^ and ([^ are 

equivalent for c = Coo • Next we obtain from p.lSp that 

iu' — ^Cooiu — w))' = u" + Ccou' u' — "fkuil — w) (4.24) 

7Coo 

and we see that (|4.22p implies (|2.ip . Conversely, from (|2.ip and (I4.24p we deduce 
that 

(u' - 7Coo(u - w))' = (u' - 7Coo(w - w)), 

7Coo 

whose solutions are 

(m' — 7Coo(w — w)) (a;) = ae~ "f"^ ^ for any a € R. 

The condition that u(x),w{x) converge exponentially to unity as a; ^ —oo implies 
that a = and (li?^ therefore holds. □ 

This reduction allows us to prove the existence of a travelling wave with speed 
Coo for (fLTjl . pr2)) by a phase-plane analysis for (|422| . (|4?23| . 

Proposition 4.7. For a/Z 7, /c > t/iere exists a monotone travelling wave 
(coo, uq, Wo) /or p.ip . (|1.2p . As x tends to infinity, uq,wo decay exponentially 
to zero with decay rate 

Aoo = Aoo(fc) := Coo7(^ - y^ + fc(l + 7)) < 0. (4.25) 

Proof. The system (li^ . has the two stationary points (0, 0) and (1, 1). We 

define the set 

J := {(a;i,a:2)^ G : < a;i < < 1} (4.26) 
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W 



1 




- stable manifold at (0, 0) 
<- direction field {u',w') 



1 



u 



Figure 1. Schematic of the phase portrait of (I4.22p . (|4.23|) . The 
dashed area represents the set J, the incoming arc the stable man- 
ifold at (0, 0) and the arrows the direction field. 

and observe that R'^ \ J is an invariant region for (|4.22p . (I4.23P (see Figure [T|). 
Next we consider the linearisation at (0, 0) of (I4.22p . (|4.23p which is given by 



for ^ : K ^ R-^. The eigenvalues of this linear system are Aoo, as defined in (|4.25p . 
and 



One checks that vi < vi holds for the components of v and deduces that the 
eigenspace corresponding to Aoo intersects the set J defined in (|4.26p . 

By the stable manifold theorem, see for example [PerOl] , Theorem 2.7, there 
exists a trajectory that, taking t as its parameter variable, converges to (0,0) as 
t —* CO and starts at t = in J, since the stable manifold is tangential to the 
eigenspace corresponding to the eigenvalue Aoo- Following this trajectory back 
with decreasing t we cannot leave J, since otherwise the trajectory would stay in 
M'^ \ J as i increases and thus could not reach any point in J at i = 0. In J, with 
decreasing the trajectory has to be monotone in both components and therefore 
has to approach the stationary point (1, 1). Thus the trajectory connects (1, 1) to 
(0,0) and satisfies (|4.22p . (|4.23p . By Lemma l46l this shows that (coo, uo, ''i'o) is a 
monotone travelling wave for p.ip . p.2p . Moreover uo(a;), w^^ix) are, for sufficiently 
large x, in the stable manifold of ()4.22p . (|4.23p and we deduce that they converge 
exponentially fast to zero, with decay rate given by Aoo- Q 



c 





An eigenvector with eigenvalue Aoo is given by 
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The existence of a travelling wave with speed Cqc implies immediately the fol- 
lowing estimate. 

Corollary 4.8. The minimal speed Cmin{k) of travelling waves satisfies 

Cmin{k) < Coo- (4.27) 

As we will see, for sufficiently large values of k the minimal speed is identical to 
the value Coo- 

4.3. Analysis of the decay- rates at -|-oo. In this section we further investigate 
the deacy of travelling waves to zero as x tends to +00. With this aim we analyze 
the functions Ck defined in Lemma 14.21 Cfe(A) is the speed of a travelling wave with 
decay rate A at +00. 

Remark 4.9. Corresponding to the reduction ()4.22|) . ()4.23p of the travelling- wave 
system (|2.2|) . we find that the equation ()4.2p for the possible decay rates 

factorises for the speed Coo = (7(1 + 7))^^^^- The value 



-— = -\/^ (4-28) 
7C00 V 7 

is for all fc > a negative root of (14. 2p with c — Coo- The decay rate of uq, wq is 
given by the other negative root of this equation, which is the value Acx) {k) defined 
in (14:2511 . 

Lemma 4.10. For k > consider the values Aiin(fc), ciin(fc) as defined in Lemma 
\4-S\ and the values Aoo(fc), A^ as given in (I4.25p . ()4.28p . Then there exists a unique 
ko > 0, which is is given explicitely by ko — {1 + 27)7^^, such that 

Aii„(fco) = a:,. (4.29) 

Moreover, for < fci < /cq < ^2 

Aoo(A:i) > Aiin(A:i) > A:,, (4.30) 

Aoo(fco) - Aiin(fco) = a:,, (4.31) 

Aoo(A:2) < Xunik2) < X*oo (4.32) 

hold, see Figure\^ 

Proof. We have proved in Lemma 14.21 that the functions attain their minimum 
Clin (A;) at a unique value Aiin(A:), 

Clin(fc) = Ck{\v.n{k)), (4.33) 

Cfe(A) > ciin(fc) for all A < 0, A 7^ Aiin(A:). (4.34) 
For convenience we recall that 

3AL = -(1 + ik) + 2v/(l + 7A)^+3fc, Aii„ < 0, (4.35) 



Aiinciin - (1 + ik) - V(T+W+3fc. (4.36) 
By Proposition [1771 Remark 14.91 and the definition of c^ in Lemma W?^ 

CkiKo) = Cfc(Aoo(fc)) - Coo for all fc > (4.37) 
holds. Next we see from (|4.35l) that Aiin(fc) is strictly decreasing in fc and that 

lim Aiin(fc) ~ ^o"^ ^ ^™ Aiin(fc) = —00. 

k — ^0 3 k — >oo 
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Since < — 1 by (I4.28P there is a unique value ko such that ()4.29p holds. By 
(gSll), (j433l) and (|4:37l) we deduce that 

Clin(fco) = Cfc„(Aiin(fco)) = CkaiXlo) = Coo, (4.38) 

and by ^f^^ that 

Clin(^o) = Cfeo(Aoo(fco))- (4.39) 

By (11331), and g^Sl) this yields 

Finally one derives from (|4.29p . (|4.36p and (|4.38p that ko is given by 

_ 1 + 27 

— 7^ — . 



To prove the inequalities ()4.30p . (|4.32p we first observe that ()4.37p implies that for 
k ^ ko there is a A between A^ and \oo{k) such that c',^{X) = 0. By ()4.10p . ()4.1ip 
we conclude that A = Aiin(fc) and deduce that Xunik) lies between A^ and Aoo(fc). 
Since Aiin(fc) is monotonically decreasing in k we obtain from l|4.29p that 

Aoo(fci) > a;„, Aoo(fc2) < a;;^, 

which proves (lOO]) . ([432]) . □ 
The conclusions of Lemma 14.101 are illustrated in Figure [H 



Ck 




Aoo(fcl) 



k = ki < ko 

k = ko 

k = k2 > ko 

(Aiin(fcl),Qm(fcl) 
{>^Hn{ko) , CHn{ko) 
{Min{k2),CUn{k2) 



A 



Figure 2. Graph of the functions for different values of k and 
the decay rates Aoo(fc) of the travelling waves with speed Coo as 
found in Proposition l4.7l 



By the previous results we can now compare the decay of two different travelling 
waves as a; — *■ +cxd, similarly as in CoroUarv 14.51 for the convergence to unity as 
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Lemma 4.11. Let k > ko and assume that (c,u,w) is a monotone travelling wave 
with speed c < Coc- Then, as x +00, u{x),w{x) decay slower to zero than does 
the travelling wave uq{x),wq{x) obtained in Proposition \4. 7[ 

Proof. By ()4.10p the function Cfc is monotonically decreasing for A < Aiin(A:). From 
and ((i?77)) we therefore deduce that 

Cfc(A) > Cfe(Aoo(fc)) = Coo for all A < Xocik)- 

Since we have assumed that c < Cqo this implies that the travelling wave (c, 
decays with a rate A > Xoc(k)- On the other hand, by Proposition 14.71 the rate of 
the exponential deacy of uq^vuq as x +00 is given by Aoo(fc)- D 

4.4. The nonlinear selection regime. In this section we prove that for k > ko, 
where fco is given in Lemma I4.10[ the minimal speed of travelling waves for the 
reaction-diffusion system (jl.ip , ()1.2|) is identical to the minimal speed of travelling 
waves for the Stefan- like limit problem (I3.2p - (|3.5p . In particular there is nonlinear 
selection of the minimal speed for k > kg. This result follows from a comparison 
principle which is formulated in the next Theorem. In general, invariant region 
arguments do not apply for elliptic systems, but here a shift parameter is chosen 
to play the role of the time parameter in the proof of comparison principles for 
paraboHc systems. 

Theorem 4.12. Let (ci,ui,wi) and {c2,U2,W2) he two monotone travelling waves 
and assume that ci < 02- Let Ai,A2 denote the decay rates at +00 of ui,Wi and 
U2,W2 respectively. Then Ai < X2 holds and, as x tends to infinity, ui{x),Wi{x) 
cannot converge exponentially slower to zero than U2{x),W2{x) do. 

In particular, a travelling wave has minimal speed if and only if its decay rate at 
+00 is the minimal one among all travelling waves. 

Proof. Let us assume that Ai > A2 holds. Since ci < C2 the travelling wave (ui, wi) 
converges by Corollarv 14.51 faster to (1,1) as x-^ —00 than does (w2, W2). In par- 
ticular 

(ui,wi)(x) > {u2,u)2){x) for sufficiently small X < 0. 

Since we have assumed that Ai > A2, the decay of {uijWi) at +00 is slower than 
the decay of (1*2, W2) and we deduce from Lemma |4 . II that 

{ui,wi){x) > {u2,W2){x) for sufficiently large a; > 0. 
This implies that there is a shift xq, such that for (ui, wi) = (ui, wi){. + xq) 

wi > U2, wi > W2 on R 
holds and such that there exists a xi e R with 

ui{xi) = U2{xi) or wi{xi) = W2{xi). 
From equations (|2.1[) . (|2.2p we obtain for U := U2 — ui,W :^ W2 — Wi 

= U" + cU' - {ci -C2)u'2-U + W '~-/k{l-W2)U + -/kuiW, (4.40) 
= cW' - (ci - C2)w2 + fc(l - W2)U - kiiiW. (4.41) 
Assume that ui{xi) = U2{xi), which gives 

U{xi) — maxU{x) = 
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and U'ixi) = 0,U"{xi) < 0. Then (144011 yields 

< - {Ci - C2)u2ixi) < 0, 

which is a contradiction. 
If wi{xi) — W2{xi) then 

W{x2) = maxW{x) = 

X 

imphes W'{x2) = and (I4.4ip yields 

< -{ci-C2)w[{x2) < 0, 

which is also a contradiction. Thus we deduce that Ai < A2. 

The final conclusion of the Theorem follows now by a contradiction argument. 

□ 



The comparison principle Theorem 14. 121 and Lemma 14.111 imply that for k > ko 
no monotone travelling wave exists with lower speed than Coo . 

Corollary 4.13. For k > ko the minimal speed of travelling waves is given by 

Cinin(^) — ^00 • 

Proof. Assume there is a monotone travelling wave (c,u,w) with c < Coo and let 
{coo,uo,wo) be the monotone travelling wave with speed Coo which we have found 
in Proposition l4.7l By Lemma |4?Tl] the functions {u, w) decay slower to zero at +00 
than Uqo, Woo do, which is a contradiction to Theorem l4.12l Therefore Cniin(fc) > 0^0 
holds and, recalling (|4.27p . the conclusion follows. □ 



5. Travelling waves in the fast-degradation-rate limit 

We complete our investigations by proving that the travelling waves of (jl.ip . 
(|1.2p are, for large values of fc, close to a travelling waves of the limit problem 

(El)- din). 

Proposition 5.1. Let {c,Uk,Wk), fc G N . be a sequence of monotone travelling 
waves for (jl.ip , p.2p with speed c > Cqo and 

Ufe(O) = i for all k>0. (5.1) 
Then, as k tends to infinity, 

Uk ^ U uniformly, (5-2) 

Uk ^ U in C°f (M), for all < (3 < 1, (5.3) 

Wk W m ifo^(K), for alll<p<oo, (5.4) 

where U, W is the unique travelling-wave solution of the limit problem with speed c 
and U{0) ^ 1/2. 

Proof We recah that u'^,w'^.<0, that < Ufc,u;fc < 1 and that by ((2Tf|) . ((2?20| 

-u'k < 1^ (5.5) 



TRAVELLING WAVES IN A TISSUE-DEGRADATION MODEL 



21 



holds uniformly in fc, where fi was defined in (I2.15|) . This yields the existence of 
a subsequence /ci — > oo — > oo) and monotone decreasing functions u,w with 
< u,w < 1 and u E hI^^ (M.) such that 

Mfc, in (M) for all < /3 < 1, (5.6) 

<^ ^ w' in £f„^(M) for aU 1 < p < oo, (5.7) 

Mfc. — > w uniformly in M for all 1 < p < cx), (5.8) 

Wfc^->w inL{;^(R) for aU 1 <p < oo,, (5.9) 

Integrating equation (|2.2p over R, we obtain 

Uk,{l-Wk,) = f, (5.10) 
and by Fatou's Lemma we see that 

/oo />OC 
u{l — w) < liminf / Uk-{l~Wk-) — 0, 
-oo J-oo 

which implies that 

u{l — w) = almost everywhere in R. (5.11) 

By dOl), (EH) we obtain that u(0) = 1/2 and u{x) > 1/2 for a; < 0. From ((5?TT1) 
we deduce that 

w{x) — 1 for X < 0. 
The equations ()2.1|) and ()5.5|) . (I5.10|) yield the estimate 



(1-UfeJ = -4^(0)-cufc,(0) + c+ / (l-wfcj+ / 7fcufc^(l-WfcJ 

D J —OO J — OO 

< C{c,j) (5.12) 

and we obtain that (1 — m) G ((— oo, 0)) and u{x) ^ 1 as x ^ — oo. By p.ip and 
(lESl), (ICTOl) 

{uk,-Wk,) = -w'j,^ (0) - - - 7fci / Mfc,(l-WfcJ > -C(7,c) (5.13) 

^ Jo 

holds and we further deduce from (|5.10p that 

{uk, - -zWk,) < / Ufc,(l-u;fe,) < — , (5.14) 
^ Jo 

which gives, substracting (|5.13p . 

oo 

Wk, < C(7,c) (5.15) 
for fc > 1. Further we find from (|5.14p that 

uk, < C{j,c) (5.16) 



holds. By Fatou's Lemma this implies that u^w £ i^((0,oo)) and u(x),w(x) — > 
as a; — > oo. By ()5.6p and since the limits as x ^ ±oo of u^^, U coincide, we deduce 
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([5:21) from |Die77[ Lemma 2.4]. 

The equations ()2.ip . (12. 2p yield that Uki,Wki satisfy 

= / (life. + c(ufe, + jWki))v' + [uki - Wki)v 



for all 77 e C^(R) and, according to (|5. 6^ - 1)5. 9p . we can pass in this equation to the 
limit i 00. This yields 

= {u' + c{u + ^w))'q' + {u-w)r]. (5.17) 



Since u, w are monotone decreasing from unity to zero and satisfy (|5.1ip . we deduce 
that there is a € M such that 

u{x) = for X > xq, 

u{x) > 0, w{x) = 1 for a; < xq. 

Therefore ([CT7|) yields that 

= u" + cm' — u + I for X < Xq, 

= jcw' + w for X > Xq, 

and that the jump condition 

7c(l — lim w{x)) — — lim v! {x) 

x\,XQ x'\xa 

has to be satisfied. This shows that u, w is a travelling wave with speed c of the 
limit system jS^), □ 

6. Conclusions 

We conclude our investigations with a brief summary and discussion of our re- 
sults. 

6.1. Slow and fast decay at the threshold c = Coo- Our results on the existence 
of travelling waves for the system (|l.ip , p.2p and the selection mechanism of the 
minimal speed are summarized in Figure [31 As in the preceding figure we have 
plotted the functions Ck which give for a travelling wave with speed c the possible 
rates A of the exponential decay to zero at 2: — > +00. By the circles, squares and 
diamonds in Figure [3] we have indicated the decay rates which in fact are realized 
by a travelling wave: presuming that there is linear selection for k < ko travelling 
waves exists for all speeds c > cun. By Theorem 14.121 the decay rates of these 
travelling waves corresponds to values on the increasing branch of (A, Cfc(A)). This 
behaviour changes for k > kg: the decay rate of uo,wo is on the decreasing branch 
of the solution curve. By Corollary 14.131 Coo is the minimal speed and the decay 
rates of travelling waves with larger speeds are on the increasing branch. 

To give an explanation of what happens in the nonlinear selection regime A: > fcg if 
the speed c falls below the minimal speed c„iin(fc) = Coo we consider the linearization 
of pTTjl . at {u',u,w) = (0,0,0), which was given in (g^)- For aU speeds c 

in a neighbourhood of Coo the linearized system has two negative eigenvalues; the 
stable manifold for (jl.ip . (|1.2p at {u',u,w) = (0,0,0) is two-dimensional. One 
checks that the eigenspace corresponding to a negative eigenvalue intersects with 
the set {z = (01, Z2, Z3) : zi < 0, Z2 > 0, Z3 > 0}. A monotone travelling wave exists 
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■oo 



k = k2 > ko 
- realized (fe = 
realized (fe = 
o realized {k = 



h < ko) 
ko) 

k2 > ko) 



Aoo(fc2) 



'OO 



Aoo(^l) 



A 



Figure 3. Graph of the functions for different values of k and 
the points on the graph that are realized by a travelling wave. 

for c > Co and approaches (0, 0, 0) tangentially to the eigenspace corresponding 
to the larger negative eigenvalue ('slow decay'), see Figure [S] We expect also for 
c < Coo an orbit {u',u,w) connecting (0,1,1) with (0,0,0). Such an orbit will 
also converge to (0,0,0) tangentially to the eigenspace corresponding to the larger 
negative eigenvalue but comes from the 'wrong' side, taking negative values for 
w, w. For the threshold c — c^c there still exists a monotone travelling wave; this 
travelling wave approaches (0,0,0) tangentially to the eigenspace of the smaller 
negative eigenvalue ('fast decay'). 

6.2. Discussion. Travelling waves often determine the long-time behaviour of so- 
lutions for the initial- value problem for p.ip . (|1.2p . Typically, for solutions with 
sufficiently fast decaying initial data, the propagation speed of pertubations from 
the unstable equilibrium is given by the minimal speed of travelling waves. The 
proof of such a result, as well as the uniqueness of travelling waves, for the system 
(jl.ip . (|1.2p is not in the scope of the present article. Nevertheless our result that the 
travelling wave with minimal speed has the fastest decay supports that conjecture. 

For the reaction-diffusion system (jl.ip . (jl.2p in arbitrary space-dimension n our 
analysis yields a family of supersolutions: Consider for c > Coo a travelling wave 
u, w with speed c and an arbitrary real number r G R. The functions iZ, w defined 



by 



u(x,t) :— u{\x\ — ct — r) , w(x,t) :— w{\x\ — ct ~ r) 



satisfy 



df U — Au + u — w jku{l — w) 
dtw — ku{\ — w) 



(6.1) 
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Since u' is negative [u, w) is a supersolution. In order to construct a subsolution 
one has to control the dimension-depending correction term in ()6.ip . 

In view of the appHcations, the robustness of the wave-speed to changes in the 
parameter values is a valuable feature: see the explicit formula (j3.15[) for the mini- 
mal speed and Proposition l5.ll For the mathematical analysis of reaction-diffusion 
systems and the selection of the minimal speed the model that we have derived is a 
good paradigm. We prove that nonlinear selection occurs and determine explicitly 
the minimal speed. One crucial ingredient is the exact first integral obtained in 
Lemma [4.6[ which is a special property of the system (|1.2p . Other results, 

in particular the comparison principle Theorem 14.121 and the observation that the 
fastest decay is realized by a travelling wave with minimal speed, can be extended 
to general monotone systems. 
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